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Abstract We investigate the half-filled Hubbard model with spatially alter- 
nating interactions by means of the two-site dynamical mean-field theory. It is 
found that a single Mott transition occurs when two kinds of interactions are 
increased. This implies that the different interactions are essentially irrelevant 
at the critical point. The nature of the Mott states is also addressed. 



1 Introduction 

Ultracold atomic gases have attracted much interest since the successful re- 
alization of Bose-Einstein condensation in a bosonic *^Rb system [1 . Due 
to its high controllability in the interaction strength, particle number, and 
other parameters, many remarkable phenomena have been observed such as 
the Mott transitions in the bosonic and fermionic systems [H[31I3]. Recently, 
the spatial modulation of the interaction has also been realized in the ^''"'Yb 
gas system [5 , which stimulates further theoretical investigations on particle 
correlations in the ultracold atomic systems. 

One of the interesting questions is how the spatial modulation in the in- 
teractions affects ground-state properties in the fermionic optical lattice. Al- 
though the optical lattice system with uniform interactions has been stud- 
ied theoretically and experimentally, the system with alternating interactions, 
which may be one of the simplest systems, has not been discussed so far. There- 
fore, it is desired to clarify how the alternating interactions affect ground-state 
properties, in particular, the nature of the Mott transition. 



T. Saitou and A. Koga 

Department of Physics, Tokyo Institute of Technology, Tokyo 152-8551, Japan 
E-mail: koga@phys.titech.ac.jp 

A. Yamamoto 

RIKEN, Advanced Institute for Computational Science, 7-1-26, Minatojima-minami-machi, 
Chuo-ku, Kobe, Hyogo 650-0047, Japan 



2 



Takamitsu Saitou ct al. 



Motivated by background, we consider the infinite-dimensional Hubbard 
model with alternating interactions. By means of the two-site dynamical mean- 
field theory (DMFT) [6] , we discuss how particle correlations affect the metal- 
insulator transition at half-filling. The nature of the Mott phases is also dis- 
cussed. 

The paper is organized as follows. In Sec.[5J we introduce the model Hamil- 
tonian and briefly summarize our theoretical approach. We demonstrate how 
particle correlations affect ground state properties in Sec. [31 A brief summary 
is given in the last section. 



2 Model and method 

We study the Hubbard model with alternating interactions, as 

H =-t^ cLcja + i^it^^l - \ ["it + niilj , (1) 

{i,j)iy a.iSa 

where c\^{cia) creates (annihilates) a fermion at the z-th site with spin o'(=t, i 
). i is the hopping integral, and Ua is the site-dependent onsite interaction in 
the sublattice a{~ A, B). At half-filling, the systems with {Ua, Ub), {Ub, Ua), 
{—Ua, — Ub) and {—Ub, — Ua) are identified since the Hamiltonian is invariant 
under the particle-hole transformations [7] as Ci-|- — >■ Ci^ and qj, — >■ (— l)'c|, . 

To investigate the Hubbard model eq. ([T]), we make use of DMFT [8l[ni 
[T0l[m[12] . which has successfully been applied to various strongly correlated 
electron systems. In the framework of DMFT, the lattice model is mapped to 
an effective impurity model, where local particle correlations are taken into 
account precisely. The Green function for the original lattice system is then 
obtained via self-consistent equations imposed on the impurity problem. 

In DMFT, the Green function in the lattice system is given as. 



G.(k,z) 



-fJ.A-SAa{z) -£k 

-Ek Z + flB-^Baiz) 



(2) 



where fj,a and Ua{z) are the chemical potential and the self-energy for the 
a-th sublattice, and ek is the dispersion relation for the bare band. In terms 
of the density of states (DOS) pa{x), the local Green function is expressed as, 

GaAz) = r r^^^, :,Po{s)de, (3) 



r]Aa{z)riBa{z) 
where 

rjaaiz) ^ Z + Sacr{z). (4) 



In the following, we use the semicircular DOS with the half bandwidth D 

2 

tD 



Po{x) = \/^ ^ i^/D)'^ , which corresponds to the infinite-coordination Bethe 
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lattice. Then the self-consistency condition for the sublattice a is represented 
as, 

^X2 



Goaaiz) = Z + fia- — j Gaa{z), (5) 

where Qq^ is the non- interacting Green function of the effective impurity model 
for the sublattice a. 

There are various numerical methods to solve the effective impurity prob- 
lem. To discuss the Mott transitions in the Hubbard model at half-filling, we 
use here the two-site DMFT method ^ , which provides us with a transparent 
view of the phase transitions [T3J[TS]. In this scheme, we introduce a specific 
Anderson impurity model for the sublattice a, which is connected to only one 
host site, as 

= 4'^ E dLda. + E + E('^-/- + 

a a a 

+ Ua (^nfc^nfai - ^ [Ufat + nfal]^ , (6) 

where daa{faa) annihilates a fermion with spin a in the effective bath (lo- 
cahzed band) and njacr = flafaa- Note that the effective parameters in the 
impurity model such as the spectrum of host particles e^a"' , the energy level of 
localized band Sa and the hybridization Vq, should be determined selfcon- 
sistently so that the obtained results properly reproduce the original lattice 
problem. When the system is half-filling, ei'''' = e^/'' — 0. 

This simple method provides sensible results as far as low-energy properties 
are concerned. Namely, for low-energy excitations around the Fermi surface, 
the Green function for /-fermions may be approximated by a single pole as 

Go.(z) ^ ^, (7) 
z 

where the residue Waa corresponds to the quasi-particle weight for the sublat- 
tice a. Therefore, the self-consistency condition is simplified as 

Vc. = (8) 

By estimating the quasi-particle weight for each Anderson impurity model, 
we proceed to perform the DMFT iterations. In the following, we discuss the 
metal-insulator transition in the system with alternating interactions. 



3 Results 



Let us discuss the stability of a metallic state in the Hubbard model with 
spatially-modulated interactions. By making use of the two-site DMFT method. 
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we calculate the quasi-particle weight for each sublattice, which is defined as 

..= (l-^^^^)". (9) 
The obtained results for fixed ratios Ub/Ua — 1.0,2.0,-0.7 and —2.0 are 



(a)f/«=(7, {h)U, = 2U, 




Fig. 1 Solid (daslied) line represent the quasi-particle weight zj^^zg) for the A{B) sublattice 
obtained in the systems with fixed ratios Ub/Ua = 1-0 (a), 2.0 (b), -0.7 (c) and -2.0 (d). 



shown in Fig.[T] When Ub — Ua, the system is reduced to the Hubbard model 
with the homogeneous interactions. As the interactions are introduced, the 
quasi-particle weights are decreased, as shown in Fig.[T](a). Finally, the quasi- 
particle weights vanish, where the Mott transition occurs at Ua = 31? in the 
repulsive case and the pairing transition occurs at Ua = —3D in the attractive 
case. These are consistent with other numerical studies such as the numerical 
renormalization group {Ua = 2.9AD) [TBI and the exact diagonalization {Ua = 
-2.98D) \Xl. 

When the spatial dependence of interactions is introduced, the quasi-particle 
weights depend on the sublattices, as shown in Figs.[T](b), (c) and (d). Namely, 
the strong renormalization appears in the sublattice with a larger local inter- 
action. Nevertheless, the quasi-particle weights vanish simultaneously in the 
strong coupling region. This suggests that particle correlations induce a single 
phase transition. To clarify the nature of the single phase transition, we fo- 
cus on low energy properties around the Fermi level. When the heavy metallic 



Title Suppressed Due to Excessive Length 



5 



state in the vicinity of the critical point is considered, the quasi-particle weight 
in the effective Anderson model eq. (O can be expanded in Va [13: as 

In this case, the effective hybridization for the other sublattice is comparable 
since it is given as Va — jp-Va- Therefore, the fermions in both sublattices are 
simultaneously renormalized and different interactions are irrelevant around 
the critical points in this approximation. The critical interactions for the tran- 
sition are then given as 

UaUb = {SDf . (11) 
The phase diagram is shown in Fig. [2] When the interaction is weak, the 
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Fig. 2 Ground state phase diagram for the Hubbard model with spatially-modulated in- 
teractions. 



normal metallic state is realized. On the other hand, when the interactions in 
both sublattices are large, there exist four kinds of the Mott states, which are 
denoted as the Mott states (i), (ii), (iii), and (iv), as shown in Fig. [2] We find 
that the Mott state (i) [(iii)] in the repulsive (attractive) Hubbard model is 
widely stabilized away from the symmetric limit (Ua = Ub)- In addition, the 
Mott states (ii) and (iv) are also stabilized in the system where interactions 
are alternately arranged repulsive and attractive. 

To clarify the nature of these Mott states, we also calculate the double oc- 
cupancy Di{= {ni-fUii)). In the non- interacting case {Ua = Ub = 0), the quan- 
tity is quarter since the empty, single occupied and doubly occupied states are 
equally realized at each site. Introducing the interaction, the quantity for each 
sublattice varies, as shown in Fig. [3] When Ub = 2Ua, the double occupancy 
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for the sublattice B rapidly decreases since the local Coulomb interaction is 
large. Finally these quantities for both sublattices vanish simultaneously at the 
Mott transition point. On the other hand, the double occupancy approaches 
half as the attractive interaction increases. In this case, the pairing transition 
occurs at Ua = —2.12D {Ub = '^Ua), where the single occupied states are 
forbidden and the empty or doubly occupied state is realized at each site. It is 
also found that when the repulsive and attractive interactions alternate in the 
system (UaUb < 0), the Mott and pairing transitions occur simultaneously in 
the corresponding sublattices, as shown in Figs.|3](c) and (d). 

In this paper, we have clarified the nature of the Mott transitions in the 
system with alternating interactions, by restricting our discussions to the para- 
magnetic case. Although this condition may be relevant for the optical lattice 
experiments, some ordered states are naively expected at very low tempera- 
tures. Therefore, it is necessary to clarify how the magnetically ordered state or 
the superfluid state is stabilized in the system only with repulsive or attractive 
interactions [UaUb > 0). On the other hand, the Mott transition should be 
realized in the system with repulsive and attractive interactions (UaUb < 0) 
since any ordered states are hard to be stabilized. In the case, the hole doping 
effect is one of the important problems since the local particle density varies 
in the optical lattice system with a confining potential. Furthermore, it is ex- 
pected that particle correlations induce the commensurability in the doped 
system |19ll20j . which is now under consideration. 
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4 Summary 

We have investigated the half-filled Hubbard model with spatially-modulated 

interactions by means of the two-site dynamical mean-field theory. It has been 
clarified that a single Mott transition occurs when two kinds of the interactions 
are increased. We have determined the ground state phase diagram and have 
discussed the nature of these Mott states. 
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